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ABSTRACT 


A portable data acquisition system was developed for 
experimental investigations of aerodynamic models. The system 
records the test data in digital form on a standard cassette tape 
recorder, which is later played back for analysis by computer. The 
early conversion of data from analogue to digital form maintained 
overall accuracies to + 1-percent. The system is notable in that it 
eliminates the need for any manual data handling, and is not restricted 
to the present application. 

The data acquisition system was used in wind tunnel and 
free flight testing of low Reynolds number airfoils. Some suggestions 
for improvements to the system arise from a discussion of these 


tests, and comparison to theoretical results. 
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CHAPTER I 


INTRODUCTION 


A large quantity of data which is mainly repetitive in 
nature must be collected and analysed when testing aerodynamic models. 
Performing this function manually is extremely time consuming and 
may easily result in errors. In the special case of free flight 
testing, it is extremely difficult to collect this data manually. A 
system which automatically collects and records series of data at a 
high speed becomes necessary. If, in addition, the system can be 
interfaced with a computer for data analysis it becomes even more 
valuable. These were the incentives for the construction of the 
system described in this paper. 

Free flight testing has been carried out with models 
mounted on trucks and airplanes (1). The advantages of this method 
of testing are that large scale models can be used, that wind tunnel 
constraints are avoided, and that free air turbulence is simulated 
exactly. The most serious disadvantages are the problems of operating 
equipment in the field and the great dependence on weather. The use 
of a data acquisition system such as is described in this paper best 
alleviates the first difficulty; the second one cannot be helped. 

There are many data acquisition systems available 
commercially. However, several restrictions and limitations placed 
on the system due to its application in free flight testing made it 


desirable to construct a system suited specifically for-these tests. 
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CHAPTER II 


DATA RETRIEVAL 


For aerodynamic testing, the minimum set of experimental 
parameters which must be measured are airspeed (or dynamic pressure), 
static pressure, angle of attack, local static pressure at points on 
the model surface, temperature, and velocity at points in the wake of 
the model. Many measurements of local static pressure and wake 
velocity must be made at each test configuration to give the pressure 
distribution and wake velocity profile. From this data, values of 
lift, moment and drag can be found, and a complete performance 
analysis can be made. 

The system developed to perform these measurements is shown 
schematically in Figure 1. The system operates under the control of 
the central electronics module. Due to the complexity of the systen, 


the major components will be discussed separately. 


2.1 TRANSDUCERS 

The transducers used for measurement of the parameters men- 
tioned above were standard models. Their general features will be 
included here, as they affect other components in the system. 

The dynamic and static pressures were measured with Validyne 
capacitance-type pressure transducers. The dynamic pressure was 
measured directly, while the static pressure, in free flight tests, 
was measured relative to a sample of air contained in an insulated 


and sealed vacuum bottle. 
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Apart from physical limitations on the size and weight of 
the system, and the necessity that it be operated under battery 
power, it was of prime importance that the system be easily operated 
with a minimum of supervision by the pilot of the carrier. It was 
necessary, therefore, to provide some kind of automatic control on 
the system. This would be necessary even if a commercial system were 
used. However, the nature of the data was such that discrete 
monitoring of several data channels (rather than continuous monitoring) 
was necessary. This sequencing function could also be built directly 
into the system controls, rather than modifying a commercial system. 
The construction of a system also ensured that the system could be 
easily interfaced with the data analysis computers. The system 
could also be constructed so as to provide digital signals recorded 
on tape which could be easily read directly into the computer. 

A commercial data system, outside of limitations on space 
and weight, would require extensive hardware modifications and 
additions for its use in the present application. Therefore, for 
reasons of uniqueness of application, assured compatability with 
other hardware, and the availability of extensive computer resources 
for analysis, the total design and construction of a data acquisition 
system was undertaken. 

The two chief goals of the designed system were that the 
process of data recording and analysis be automated as much as 
possible, hence no manual processing of data be required, and that 
the system be portable, with the application to free flight testing 


in mind, 
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The designed system is an extension of an idea presented by 
Whitfield (2), in which a much simpler data recording system used in 
the performance testing of gliders was described. A fundamental 
difference of the present system compared to Whitfield's is the use 
of standard transducers which produce a simple analog DC voltage 
output signal. Since the central electronics accepts a standard 
voltage input (usually 0 to 5 volts) any transducer having similar 
output characteristics can be used. This greatly increases the 
flexibility of the system making it useful for a wide range of 
applications. For example, studies in vibration, fluid mechanics and 
heat transfer, to name a few, could all be carried out using this 
data acquisition system. 

The system itself operates in two parts: data recording and 
data analysis. The former is dependent largely on the hardware 
system and is described in Chapter II. The data analysis phase is 
most dependent on the software package developed for the large ground 
based computer. It is discussed in Chapter III. Following the 
system description, two applications are presented in Chapter IV, along 
with discussion of results and some suggested improvements to the 


system. 
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Figure l. Data Acquisition System Schematic 
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Temperature was measured using a TSI Telethermometer trans- 
ducer. 

A scanivalve was used for measurement of local static 
pressure at points on the model surface. This instrument is a 
combination transducer and solenoid drive mechanism. Up to 48 
pressure tubes, coming from taps on the model, can be connected to 
the ports on the head of the scanivalve. Each port is connected in 
turn, by means of a channel in the scanivalve rotor, to the transducer 
cavity. The rotor is moved from port to port by the solenoid drive 
mechanism. After a complete revolution of the rotor, a limit switch 
indicates that all the ports have been monitored. This instrument 
was ideal for this application because of its compactness and simplicity 
of operation. 

The angle of attack was measured, in the free flight tests, 
by a yaw probe mounted on the model support structure. The scanivalve 
was used to monitor the pressures coming from the two sides of the yaw 
probe. 

Two types of transducers were used for the measurement of 
velocity at points in the wake. The first was a constant temperature 
hot film anemometer. This instrument was used in hopes of getting 
information on the free stream turbulence level, which is of some 
importance to low Reynolds number airfoils. However, extreme 
difficulty was experienced in using this instrument, both in operation 
and analysis of data. Quantitative information on turbulence levels 
could not be obtained directly with the type of recorder used. An 
FM recorder is necessary for this. Also, since the anemometer had to 


be calibrated for both temperature and density variations, the analysis 
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was quite difficult. Following these problems, a standard pressure 
transducer and Pitot probe were substituted to measure the dynamic 
head in the wake directly. This greatly increased the reliability 

of the data and simplified the analysis. The probes for these 
transducers were mounted on a stem behind the model. The stem could 
be traversed through the wake by a geared variable speed motor drive. 
The position of the probe in the wake was monitored by a 10 turn 
potentiometer coupled directly to the motor shaft. The potentiometer 
was supplied with a constant voltage from the central electronics, 


and the output voltage was a linear function of the probe position. 


2.2 ELECTRONIC DATA HANDLING AND CONTROL 

The heart of the data acquisition system is the data handling 
and control system. The function of this component is to sequentially 
monitor the transducer signals and produce the digital data signals 
sent to the recorder. The control circuitry also operated the 
scanivalve drive mechanism and responded to the start and stop signals 
from the operator control and the cycle limit switches. A block 
diagram of the central electronics is shown in Figure 2 and circuit 
diagrams are in the appendix. 

The operation of the electronics in the treatment of data 
was as follows: A multiplexer circuit sequentially gated each of the 
transducer signals to a voltage controlled oscillator. The transducer 
signals were amplified and bias voltage applied where needed to 
operate the oscillator in a higher frequency range. Strict timing 
for the multiplexing process was achieved by using an electronic clock 


based on a very high frequency (262.144 KHz.) crystal oscillator. 
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Each transducer signal was gated through to the voltage controlled 
oscillator for a period of 289.05 milliseconds. The multiplexer then 
closed off the transducer signal and, after a short delay of 39.06 
milliseconds, connected the next transducer in sequence through to the 
oscillator. The signal gap was used to separate data points as will 
be explained in section 2.5. Each transducer signal was monitored 

in turn and the multiplexing process repeated until all the data was 
recorded. 

The voltage controlled oscillator produces an output 
frequency proportional to the input voltage. Input voltages from 0 to 
10 volts produced output frequencies from 0 to 10 KHz. The oscillator 
was operated in the range of 300 to 3,000 Hz. During the gap between 
transducer signals, output from the oscillator was blanked to leave 
a clean break between data points on tape. The output of the oscillator 
was then amplified and sent directly to the recording head of the 
data recorder. 

The resulting record of data on tape consisted of a series 
of data bursts of varying frequencies separated by short gaps, in 
repeated groups of four or five. The number of bursts in each group 
was determined by the mode of operation of the system. In the wake 
traverse mode, five transducers were active. In mode two, the model 
pressure distribution was scanned and four transducers were used. 

The number of cycles in each burst was therefore a digital 
representation of the voltage output of the transducers and hence of 
the physical quantity being measured. Tape speed on recording or 


playback was immaterial, since each data burst had been written on 
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tape for a precisely controlled time period. That is, neither amp- 
litude nor frequency of the signal as it appeared on tape mattered, 

as it was the number of cycles in the burst which carried information. 
The means by which these data bursts were read back and reduced to 
useful data is discussed in later sections. 

When operating in the second mode, the scanivalve drive 
mechanism was controlled by the electronics package. Immediately after 
the scanivalve transducer signal was monitored and recorded, the 
scanivalve was stepped to the next port. The transducer was not 
monitored again for slightly more than a second, giving ample time for 
equilibrium to be reached at the new pressure. 

The completion of data monitoring for a test condition was 
indicated to the central electronics by limit switches installed in 
the scanivalve drive and the traverse motor drive mechanisms. When 
these limits were met, the multiplexing and recording functions were 
terminated and an indicator lamp on the operator's control box 
signalled the end of the cycle. 

The overall control of the system was through a hand-held 
box connected by multiconductor cable to the central electronics. 
This box held all the necessary controls for power, mode selection, 
traverse direction and speed, cycle start switches, and the end of 
cycle indicator lamp. The operator could therefore move around while 
data recording was in progress and remain in full control of the 
system. 

The electronics remained idle while the next test condition 
was set up and until the signal to start recording was received from 


the control box. 
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2.3 DATA RECORDING 

The output signal from the central electronics could be sent 
directly to a real-time computer to give an immediate read-out of 
data. However, to give portability to the system, a medium for 
temporary data storage was required. Ideally, this recorder would 
have written data directly in a suitable code to computer-compatible 
magnetic tape. This was unacceptable, however, due to the size of 
such a recorder and the required power loads. A small self-powered 
cassette tape recorder was selected as the temporary storage device. 
This necessitated a further step in data processing -- the transfer 
of data from cassette to coded magnetic tape on a small ground based 
computer. This transfer phase will be discussed in section 2.5. 

It was mentioned previously that the data output signal was 
sent directly to the recording heads of the cassette recorder. The 
recorder input signal did not require amplification and the internal 
circuitry of the recorder would have added noise and distortion to the 
signal. The recorder's playback circuits were used, however, with no 
ditticulty. 

The recorder used was a two channel model. The first 
channel was used exclusively for the test data. The second channel 
was used for a voice commentary on current test conditions, recording 
start and other information to identify the data files. Originally, 
it was planned that the second channel would also record the infor- 
mation on turbulence from the anemometer. This was impractical due 
to an automatic level control in the internal circuitry of the 


recorder. 
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The recorder was mounted so as to be easily removable from 
the instrument box. The use of cassettes for data storage was 
advantageous in the ease with which they could be manipulated and 


stored. 


2.4 POWER SUPPLY 

Power was supplied to all components in the data system from 
a set of 12 volt batteries. Two motorcycle batteries and one automobile 
battery were used. The former were connected in series to provide 12 
and 24 volts to components requiring DC power. The larger automobile 
battery was used to drive an inverter which produced 110 volts AC for 
the other components. When operating in the lab, all components could 
be operated under external power to conserve the charge in the batteries. 
The traverse motor was powered separately by a 9 volt dry cell battery. 

Under internal (battery) power, the data acquisition system 
could be run continuously for two to three hours, depending on initial 
charge in the batteries. After this time, or if the batteries were 
low, performance of the system became erratic. This was usually 
indicated first by the scanivalve solenoid drive failing to step to 
the next position, causing the read cycle to take inordinately long 
to complete. However, power failure proved to be of little problem 


as all necessary data could be recorded in much less than two hours. 


2.5 DATA TRANSFER TO CODED TAPE 
With the raw test data recorded on the cassette, the last 
step in the retrieval process was to transfer this data to computer- 


compatible tape in a usable form. This was done on a small real-time 
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computer, a PDP 11/20. A functional schematic of the transfer process 
is shown in Figure 3. 

The cassette recorder was interfaced with the computer 
through a simple signal conditioner. Each pulse from the cassette 
activated a trigger in the conditioner, sending a more uniform and 
sharply defined pulse to the computer. The threshold for trigger 
activation could be adjusted so as to eliminate noise induced trigger- 
ing. This was monitored by visual inspection of the output signal 
displayed on an oscilloscope. As each pulse reached the computer, a 
count register was incremented, so that the number of cycles in the 
burst would be obtained. 

The signal conditioner module also sensed the gap between 
data bursts. When no signal arrived within 30 milliseconds (set by 
operator), an interrupt was sent to the computer. This caused the 
current total in the count register to be stored in core memory, and 
the register to be reset to zero for the next burst. 

When a certain number of count totals, specified by the 
operator as the block length, usually 40, were stored in core memory, 
the entire block was written on nine track magnetic tape in a code and 
format compatible with the large data processing computer (IBM 360/67). 

When the data set for a test condition had been completely 
read into the computer, the counting function was terminated and all 
remaining totals in memory were written to tape. This termination 
occurred on the generation of an interrupt from the teletype console 
by the operator. The resulting file on magnetic tape consisted of 


several blocks of numbers representing the number of cycles in each 


= IGT = 


wander sehaunas its a9 stannedyatniuhd 9p0 rN 

: ° 
; un -_ “E sniegt nu muito at 3 
a oe ey, : 


as20gmop ‘as Adiew boselsotis soy xeirrose7 at TereaD, aE a 


a 


 go9sdans “els wort setin ion) =. sanuts bsaos fongie aiqmie « vigeonts 
i 7 7 - 


: 7 i. : 7 7 
bie urwlim ove s wntbas .yemotsisous of; of 19a8RET2 & hetavise 


tagatay sod Blodesaa sdT .asteqnos sis ot satug peaties ¢lerade: 


7 


~Tauga18 anda sefin Sitpinils G3.an oe bejeuihe o¢ Wiuel uO iaawisoe | 
7 a 


‘ 7 7 
feaeikn tuqie silt Yo nolsaaqanl letuly yt Gexci tom aav awit sacl 
7 Sp 
& ,Jesuqgoy sii bedess1 veluq done 24) = .aqoseoitiveas me te heyniqeth 
or 


“aia nt aelseo Ap sedivin siz zadd os yhyvosmeton! any vetedyst 3969: 


1 ioe © 
peghazde of bine dered 


—~ oO 


Gsev3es Gee aia boende chic slubom tenshotdion fapgte mi? ; 
vd $68) ahidosetilia OC atitiw Goelaus [aeQce un oe .esered ejab 
wae 


Sng beaiay oii/l deditgmys eis 6! her aby 1quasszat fie , (XOoMLSqo. 


(oaee 
bor ,viomem 2100 al barote 4g 63 teJelgey Innes at! “i latoy JasTIH 7 
. > 
5 (Jazud axon add wo? etes 03) i6uan ed oF t99eiges ord a. 
. ' wen 
: 7 ; 7 } : <a 
off Yd hetitusqe ,@leyod Joos Io Sede alhicso #& asi 
: ~~ g 
| a ea son a? hesate srow WA vilevan svidais fonts i oa2 ge roratege 
| a 
- bas aboa @ Al sqéa Sissdgem Anny} Sito) 1 nee yi7 sf SAW dsofd ark de - oat 


(valor wi) Joti. Bol vasow — agxel an ftw tastes see 
, ee 7 
— vissalyae med bad poltibtoes: J#s2 # 107 a deh 4: iW ead ' 


panties? apw nolsonid dks ats ic aeons 


O ee. 


_ aa : : : . 
- | anue nal ww éoqed o3 idan aie ad sei acini aig eta 
at 2, ae ie mera agirtasn! a8 0: . 


14 
CASSETTE 
RECORDER 


TRIGGER 
THRESHOLD 
ADJUST 


| SIGNAL CONDITIONING 
AND GAP DETECTION 
(HARDWARE) 


PULSE 
TRIGGER 


| 

| 

| 

| 

| 

| 

| 
ei COUNTING MODE 

OPERATION OF 

PDP 11/20 
(PROGRAMMED) 


INCREMENT 
COUNT 
REGISTER 


STORE COUNT; 
RESET COUNT 
REGISTER 


BLOCK 
LENGTH 
Le 


WRITE BLOCK 
TO TAPE 


Figure 3a. Signal Conditioning and PDP 11/20 Counter (Schematic) 
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PDP 11/20 Count Initiation and Interrupt Servicing 
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burst of the cassette record read. These numbers formed a data 
vector, with no distinction between each number to indicate its 
transducer origin. The job of sorting these numbers into separate 
vectors was done in the programs described in Chapter III. 

After the terminate interrupt had been processed, the 
computer prompted on the teletype for the desired blocksize and file 
label for the next data set. When the counting function for the next 
data set was initiated through the teletype, this file label was 
written on the magnetic tape as the first block in the new data file. 
The label could be recovered on the large computer to give details of 
the test conditions. 

The program for the PDP 11/20 is written in Assembler and 
is included in the appendix. This program resided on a magnetic tape 
in core image form and was loaded into the computer before each data 
transfer session. 

With the digital representation of data on tape, the next 
phase in the data analysis could be carried out on the larger data 
processing computer. This was largely a software dependent function 


and is described in the next chapter. 


2.6 CALIBRATION AND ACCURACY 

Before any testing was carried out, the system had to be 
calibrated. This was done by setting up known physical inputs to the 
transducers and operating the system as if test data were being 
recorded. Using a small auxiliary program which simply read the 
coded magnetic tape, a list of counts for each transducer for each 


input condition could be obtained. These count totals were used to 
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determine the calibration relations, samples of which are shown in 
Figures 4 through 9. The use of the derived calibration relations is 
discussed in more detail in the next chapter, but it should be mentioned 
here that the system was entirely recalibrated before every test 

series. This procedure required approximately half an hour. 

The accuracy of the system can be described in three ways: 
precision, sensitivity, and repeatability. The precision of the data 
system refers to how accurately the computed value of the physical 
quantity compares with the actual value. Sensitivity relates to the 
smallest perceptable change in the computed value as a percentage of 
the full range of values. Repeatability refers to the degree with 
which the same measured value results from repetitive measurements of 
the same physical quantity. 

The precision of the data system depends almost entirely on 
the calibration procedure and the resulting calibration relation. 
However, two factors other than the exactness of the calibration 
affect the precision. The first is the fact that the analog signals 
from the transducers were, effectively, integrated over the gate- 
open time. Thus, small variations in the physical quantity were 
averaged over that time period. Secondly, the count total arrived 
at by the small computer was accurate to only plus or minus one count. 
In a total range of five hundred counts (an average value) this 
represents a system induced error of + 0.2 percent, which is the lower 
bound on the precision. The upper bound on the precision depended on 
the care and accuracy with which the calibration was accomplished. 


This is estimated to be on the order of 1 percent of full scale. 
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In this situation, sensitivity and repeatability are related. 
It was found that repeated measurements of the same physical quantity 
resulted in count totals that differed by, at most, + 2 fromthe 
average. Thus, the smallest significant perceptable change is pro- 
duced by a change of three in the count total. This represents a 
sensitivity of + 0.6 percent of the full scale of five hundred. Since 
the sensitivities of each transducer were much lower than this, the 
+ 0.6 percent of full scale repeatability was the limitation, also, to 
the system sensitivity. 

The comparison of the sensitivity and precision shows that 
even if great care was used in obtaining the calibration relations, 
the overall sensitivity of the system limits accuracies to + 0.6 
percent of full scale values. However, in view of the fact that final 
results depended on integrations of these values (for example, the 
lift coefficient derived from the pressure distribution), the accuracy 
of the final results were estimated to be on the order of + 1 percent. 


This is quite good accuracy for experimental testing. 
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CHAPTER III 


DATA ANALYSIS 


Starting with the magnetic tape records of count totals, 
the task of data analysis breaks down into five main parts. The 
first is to sort the string of numbers on tape into vectors, one for 
each physical quantity being measured. Next, using suitable calibration 
relations, the actual values of the physical quantities can be cal- 
culated from the count totals. With these basic values, other physical 
quantities can be derived and secondary parameters calculated. Fourth, 
the results are made available for any further analysis desired by 
storing them on disk. Lastly, final analyses such as lift and drag 
calculations and plotting can be carried out. This last step was 
carried out by separate programs described below. 

The first four steps are included in both the programs 
SCANREAD and TRAVREAD which are listed in the appendix. A flowchart 
for these programs is shown in Figure 10. The complete specifications 
for input format, peripheral device assignment, and execution commands 
under MTS (Michigan Terminal System) control are included in the 
appendix listings. It was necessary to use two programs for this 
first phase of data reduction because of the differing formats of data 
for the two operating modes of the recording system. That is, mode 
one (traverse) records five physical quantities, while mode two (scan) 


records only four. 


Provision was made in these two programs for calibration 
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constants to be included with each set of data rather than being built 
into the program. This was done by placing the appropriate factors in 
a disk file and reading them into the program during execution. 
Alternatively, the values could be entered through the terminal, if 
that was the *SOURCE* default. The derivation of the calibration 
relations was described in section 2.6. All the relations (with the 
exception of the hot film anemometer calibration) were linear. The 
anemometer required a two stage reduction, both of which were non- 
linear. It was found that the linear calibrations were most affected 
by a zero shift, while the proportionality constant stayed the same. 
Thus, by obtaining the count values at zero conditions only, a rough 
but reasonably accurate calibration could be obtained. 

Once the primary physical values had been determined using 
the calibration relations, they could be used in combinations to 
derive the secondary values such as airspeed, density, altitude and 
pressure coefficients. 

The fourth step mentioned above, and the last operation in 
SCANREAD and TRAVREAD, was to output these results to the printer and 
to transfer them to disk for further analysis. Accompanying this data 
to disk was a list of parameters describing the test condition, for 
example, airfoil name, angle of attack, location of static pressure 
taps on the model, and flap configuration. 

Several programs were used for final analysis of the data 
written to disk by SCANREAD and TRAVREAD. The program SCANPLOT uses 
output data from SCANREAD to plot the pressure distribution on the 


model. The program SCANLIFT, using the same data as SCANPLOT, calculates 
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the lift and moment coefficients. Each of these were set up to 

process several complete sets of data in a single run. The data set 
was a compounded list obtained by the execution of SCANREAD on several 
raw data files. The program TRAVDRAG reads a compounded data list, 
produced by several runs of TRAVREAD, to calculate the drag coefficient. 
Flowcharts for these analysis programs are shown in Figures 11 through 
P35 

All the above-mentioned programs were stored in object form 
on disk. Once a data tape had been produced, the task of analysis 
reduced to issuing a set of instructions for the computer to run a 
certain program with a certain data file. The programs are entirely 
compatable with each other regarding internal data exchange formats, 
and instruction queuing. In only a few minutes, a set of instructions 
could be assembled which would produce all the desired final results, 
including complete listing a raw and reduced data, pressure distribu- 
tion plots, and lift, drag and moment coefficients. 

The programs as written rely heavily on the very efficient 
and extensive peripheral device management programs resident in MTS. 
(In running SCANREAD, seven devices were used for input/output 
operations.) Also, extensive use was made of system and mathematical 
library subroutines, some of which may not be available on all 
installations. The resident subroutines used are listed in the 
appendix. 

This completes the description of the data acquisition system. 
The next chapter describes two applications of the system, followed by 


a discussion of results and suggested improvements. 
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CHAPTER IV 


APPLICATIONS OF THE DATA ACQUISITION SYSTEM 


4,1 TEST CONFIGURATIONS 

The data acquisition system was used in two applications with 
some minor variations in power supply and transducer arrangements. The 
first was free flight testing of an FX-61-163 profile at low Reynolds 
number. The second was wind tunnel testing of the same profile with 
a 35 percent chord slotted flap adjustable to any position. 

The free flight tests were performed on a two dimensional 
model having a span of 48 inches and a chord of 18 inches. The model 
was mounted between rectangular endplates 54 inches above the wing of 
the carrier sailplane. The angle of attack was adjustable from the 
cockpit by means of a pushrod and lever mechanism, With this model 
configuration, Reynolds numbers based on model chord of from 0.6 x 10° 
tow lex 10° could be achieved. The data acquisition system was 
mounted in the rear seat of the sailplane, with the control box in 
the cockpit accessible to the pilot. 

The normal test procedure was to tow the carrier sailplane 
to 3,000 feet above ground, during which time measurements could be 
taken. Following release, more measurements were taken as the sail- 
plane descended in free flight. An entire flight lasted from 25 to 
35 minutes allowing a complete range of angles of attack to be covered. 

A very large number of problems were faced in performing the 
flight tests. Generally, these arose from the operation of a fairly 


sophisticated electronic system far from the usual support equipment 
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available in the laboratory. Numerous problems arose in the develop- 
ment of the system itself, many of which did not appear until after 
the system was mounted in the sailplane or had been aloft for the test. 
Since the airfield was some distance from support equipment, a great 
deal of time was lost in dismantling and transporting the system 
between the airfield and the laboratory. The dependence on weather 
proved to be another great problem. 

Specifically, the greatest difficulty with the free flight 
tests was the operation of the hot film anemometer. Originally planned 
to be used to get an indication of the free stream turbulence as well 
as measuring the wake velocity, the anemometer required much more 
attention than was possible in a free flight test. Also, the output 
signal during the tests was found to be well outside the range of 
values obtained in the laboratory, upsetting the electronics consider- 
ably and destroying many sets of data. This was probably caused by 
the extreme deviations of air temperature and density in the test 
conditions from those of the calibration conditions. This posed no 
great problem, however, since the wake traverse data could be obtained 
using a Pitot probe and a standard pressure transducer. This sub- 
stitution was then made on later tests. Since it was only hoped to 
get qualitative information on free stream turbulence with the original 
system (quantitative data would have required much more equipment, and 
be even harder to operate), the loss of this information was not 
considered serious. 

Preliminary results obtained in the free flight tests pin= 


pointed some further problems. Foremost among these was the 
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determination of the angle of attack of the model. As mentioned 
previously, the angle of attack was measured with a yaw probe 
mounted on the model endplates, with pressure lines leading to the 
scanivalve. Thus, the angle of attack was monitored once in scan 
mode and not at all in traverse mode. Although the pressure dis- 
tributions measured during these tests were fairly smooth, there was 
a very large scatter in plots of lift coefficient against angle of 
attack. This could have been caused by changes in the attitude of 
the sailplane. Other than the location of the pushrod, there was no 
way of determining exactly what the angle of attack was at each point 
ince cest. 

Lacking the time to make the necessary changes to the data 
acquisition system to correct these problems, and also because of 
poor weather, it was decided to terminate the free flight tests. 

For the second application of the data system, a series of 
wind tunnel tests were set up. The model tested was again the FX-61- 
163 with 35 percent slotted flap, but with only a 12 inch chord. The 
flap was adjustable to any setting. There were nearly twice as many 
pressure taps on this model compared to the flight test model to 
obtain more definition in the pressure distribution. The data system 
was run under external power to eliminate the periodic need for 
battery removal and recharging. 

It was possible to complete a series of tests in the wind 
tunnel covering a full range of angle of attack at a given speed and 
flap setting and record all the data in less than half an hour. The 


transfer of data from cassette to magnetic tape took equally as long. 
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With all the data on tape, only several minutes were required to 
assemble a set of instructions for analysis on the large computer. 


A discussion of these results follows in the next section. 


4.2 RESULTS 

The wind tunnel tests were performed at flap settings of 
0, 10, and 20 degrees, labelled Series III, IV, and V respectively. 
(Series I and II were preliminary tests used in setting up the test 
equipment.) The Reynolds number for all tests was 0.8 x 10° based 
on chord. The final results for these tests are shown in Figures 14 
to 16. A few sample pressure distributions as they are plotted by 
the computer are shown in the appendix. Several of these plots have 
been overlaid and are shown in Figures 17 and 18. 

It should be noted that all data shown is uncorrected for 
wall constraints. The use of a semi-open top, as suggested in (4), 
makes the theoretical wall corrections negligible. More experimental 
tests must be carried out to see if that is, in fact, the case. 

The present series of data show very little scatter, both 
in the pressure distributions and in the integrated lift and drag 
results, and are reasonably close to the values reported for this air- 
foil (3). Differences are on the order of 10 percent. It was felt 
that the chief sources of error in the final results were with the 
tunnel itself, since the model was quite large. This is a serious 
problem when testing high lift airfoils, and further efforts are being 
made to develop the test facilities. These efforts are directed at 
boundary layer control on the tunnel walls, and further investigation 


of the semi-open top as mentioned above. The development of this 
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facility should make it possible to investigate the airfoil performance 


at high angles of attack, going beyond the stalling point. 


4,3 DISCUSSION OF RESULTS 

It should be noted, referring to Figure 19, that the pressure 
distributions were measured and plotted separately for the upper and 
lower surfaces and do not join. Since the pressure taps did not 
extend to the trailing edge of the main airfoil, values of pressure 
coefficient had to be extrapolated in this region for the calculation 
of lift. This was done in the program SCANLIFT, to approximately 
80 percent chord. From the appearance of the plots, this should cause 
no great problem as the pressure recovery in this region is quite 
smooth. A further point of interest is that, on the flap, the exact 
location of the stagnation point can only be inferred from the data 
given. In fact, the computer-drawn plots may be misleading in some 
cases, as the plotter only connects the input data points, and cannot 
infer the location of such strategic points as the stagnation point. 
This is a minor disadvantage of a completely automated system, that is, 
the system can handle the general case flawlessly, while errors may 
arise in treatment of some specific cases. This disadvantage can only 
be overcome by providing more data points to complete the curves (more 
taps on the model), by very sophisticated programming, or by manual 
editing of results. 

Despite the fact that the results obtained from these two 
test configurations were somewhat disappointing, the data acquisition 
system was able to show its value in this type of experimental work. 


With minor modifications to the test facilities, it should be possible 
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to accurately measure the airfoil performance characteristics. 

Through the use of the data acquisition system as originally 
conceived, and an appreciation of the difficulties involved in the 
test procedures, several possible improvements to the system arise. 
None of these are major, and do not change the central idea of the 
system. 

Firstly, it would be possible to decrease the gate-open time 
of the multiplexer and at the same time increase the voltage offset 
of the transducer signals. The former would speed up the recording 
process while the latter would increase the frequency output of the 
oscillator, thereby maintaining the accuracy. The real time computer 
used to count the cycles can operate much faster than the range where 
it was used. However, for flight testing, speeding up the recording 
of data would increase the effects of instabilities of the carrier, as 
the analog signal integration period would be much shorter. 

Secondly, the method of measuring the angle of attack in 
free flight tests must be changed. It should be monitored continually 
through the test. This involves the addition of another transducer and 
the rearrangement of sequencing in the multiplexer. Also, since 
temperature is the most stable parameter, it need be monitored only 
once for each test condition. 

Since one of the ideas of free flight testing is to simulate 
atmospheric turbulence exactly, it should be measured by placing an 
RMS meter on the anemometer output, to give a DC signal for recording. 
The anemometer must also be temperature compensated. However, in view 


of the difficulties experienced with the anemometer, this would require 


very special care. 
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The process of data transfer to magnetic tape is slow due 
to the necessary delays in setting up the cassette at the desired 
position, entering file label through the teletype and starting the 
program. This could possibly be improved to be a continuous read-out. 
However, the problems of file identification and data sorting would 
be immense. With practice, it was possible to set up the computer 
in the gaps between data sets and the cassette could be run continuously. 
Lastly, the analysis programs could be improved so as to 
decrease the amount of supervision needed. It would be possible to 
queue the computer by key words in the file label. In effect, the 
programs as written would become subroutines under a master operating 
program. It would also be possible, since the small computer was not 
operating near capacity, to do the first four stages in the data 
reduction at the same time as the raw data was being read in and counted. 
However, the PDP 11/20 does not have the sophisticated subroutine 
library and input/output management program that is available on the 
larger IBM 360/67. The program for the PDP 11/20 must also be written 


in Assembler rather than FORTRAN. 


4.4 SUGGESTED IMPROVEMENTS AND FURTHER WORK 
From the preceding discussion, the following recommendations 
can be made to improve the system performance: 
1. The transducer input sequence should include an in-sequence 
measurement of angle of attack for free flight testing. 
2. If information on turbulence is desired, an RMS voltmeter should 
be used on the anemometer output to give a DC voltage signal 


which would be monitored like any other transducer. 
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3. The possibility of continuous playback of the cassette should be 
investigated more fully. 

4, The data reduction process on the large computer could be stream- 
lined to simplify the operating instructions. 

These changes, which do not drastically alter the makeup of 
the data acquisition system, would be quite easily implemented and 
would improve the performance and usefulness of the system. 

The data acquisition system has shown its usefulness for 
this kind of repeated testing. It is therefore suggested that it be 
used to make a complete study of airfoil performance in the wind | 
tunnel under many different wall and flap configurations, and boundary 
layer control. It should be possible, then, by comparing experimental 
results with theory to develop the airfoil testing facility so as to 


provide completely dependable and accurate results. 
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CHAPTER V 


CONCLUSIONS 


The data acquisition system described in this report has 
shown itself to be a useful and efficient research tool, thereby 
achieving its objective. Its operation is extremely flexible, 
allowing fixed or mobile operation,and it is not limited to aerodynamic 
testing or to the specific transducers used here. Since no manual 
processing of large masses of data is required, either in recording, 
transfer to computer punch cards or files, or analysis, a great deal 
of time is saved and processing errors are eliminated. The system 
therefore allows much more efficient use of test facilities. The 
system also combines the use of inexpensive and low powered measure- 
ments with the use of the large data processing computer. 

The present system was estimated to have an accuracy of + 1 
percent. The conversion from analogue transducer signals to digital 
representations at an early stage eliminates many sources of error 
and assures the best possible accuracy consistent with the transducer 
calibration and stability. 

It was shown that, even with automatic data recovery, the 
difficulties of performing free flight tests are quite large and 
numerous, though not insurmountable. The wind tunnel tests indicate 
that more work must be done to accurately assess the effects of the 
tunnel walls on the results. 


Data acquisition systems such as described here will become 


49 


. i“ 7 : - i 
7 - : - 
4 5 7 ' _ —— 7 
: ij : 1 
, _ 7 ” — \ : a, ; 
| eo  Weawna 
7 - - : . 


eHOTev JIWO5 


; run | 
antl 320997 #14) ni dedits2ek metayh RoLDieTapte KTER eit 


yoogsila] ,Locd dsvenesy tnstottis bas tuted « od oa ifpaak t 
: " Z : 7 a. 
9ldiveli vismesixs el rolisisqe a9t evitoetdo et} gmrivation, 


: eet! 
ilgenyborss (03 heidtwkt go 2) 21 bre, act3areqo/Si)atem to Gears en 
leu on a0gFP .sxed heew e¥sorbensyi os) tiesqe aifa o¢ 40. guiias 


; Bly 
‘gauthvoser cl asa7ts .teituphes at nish to aseeca Sgr a 6 gnleaezozg 


j@sh i#axg 6 ,@2ayleos 1s .20!22 sn sbeso. Weneq 19e0qnos 05, tele 
’ a § 
matave  syfl beitrig}is “ts erwiys goletewrag dae be ye wt sabisSen 


_ 
9iT -estaklinc? ges7 ic seu desiotits 2mm dow wolie eiolazds 


- 
; 


~Squcaem be TaWuG wot bie svlenaqyoa Jo Aza) Sit rsfifdmar ole exaatie 

-igieqnon. guieesoosq Bish Syxel Sad Y.oaw sel tsi awn 
i * to RUMOR UB sved a4 Stag lrao eri asievn soerd sAT - : 
heiioth oF alatgle racubsnt 5 aigolsn mort solexeveos mil irsoaag 


Wie. 26 UsutveS toon e5ienbmlIs syhte ‘yizse wx Ie srutiesneaszg 


r 
7 | a 6 
easubetgay sdy g5iw to9tetanos yoawwSEss sidtadog sssd ori siti 
- ay ’ - : ; ' 
: - “pyattiriere Kins feo¥den¢ aS 
7 ray @ a - 


4 


o61 ,XtaVGgQa7 BSED ih pinvodunn Hatw qevs anil inode 2aw 2 a, 7 


7 : 2 a pa 
' bee wt hens 978 #Jn9d — 5e02 spbemptreg jem ath 7 


LPs 7 a 
- wheohiad mes fatoss buy ant oldosnuoneisaly ior Pane a or 
es: 30 


mo 
PONM ie reeks 
_ - - ; - . ’ ay , pidatet 
a - _ Pe | - ; : 


jonah enseae wWayeivoas a) ere ee teu azen 
7 ifet 1 wee 


ait 


50 


more prominent in future research. With advances in technology and 
the growing abundance of minicomputers it will become increasingly 
possible to perform experimental research quickly and accurately 
without great expense. Especially, any investigation requiring the 
repeated monitoring of several data channels should be considered as 


requiring a data acquisition system such as described in this report. 
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APPENDIX A 


ELECTRONIC CIRCUIT DIAGRAMS 
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COMPUTER PROGRAMS 
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COMPUTER PROGRAMS 


Following is a list of all system subroutines used in the 


program package listed in this appendix. These subroutines were 


available in the Scientific Subroutine Package library (*SSP MATHLIB) 


or in the Calcomp plotter library (*PLOTLIB) available through MTS. 


Complete documentation for these subroutines is available in "System/ 


360 Scientific Subroutine Package" published by IBM, in "The Michigan 


Terminal System Volume 3: System Subroutines" published by the 


University of Alberta Computing Services, and in "Calcomp Users Manual" 


also published by the University of Alberta Computing Services. 


Name 


ALI 


ATSG 
AXIS 


CONVAE* 


READ 


SYMBOL 


Function 


Aitken-Lagrange interpolation at point specified, 
using ordered vectors produced by ATSG. 

sets up ordered vectors of data points for ALI. 
plots, labels and dimensions a linear axis. 
converts characters coded on magnetic tape in ASCII 
to EBCDIC which is IBM 360/67 internal code. 

plots smooth curve through data points. 

changes plot origin. 

opens a plot file. 

integration of equally spaced function values using 
Simpson's rule. 

reads data blocks from magnetic tape. 


plots specified character string. 


* This subroutine was not available through MTS and was written 
by L. Lewis, Tech. Services, U. of A. 
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58 
PROGRAM NAME: SCANREAD 


OBJECT NAME? SCANREADOB 


THIS PROGRAM OBTAINS DATA FROM MAGNETIC TAPE RECORDS 
FOR SCAN MODE OPERATION AND REDUCES THE RAW COUNTS TO 
PHYSICAL QUANTITIES. THIS VIS TE OLLOWED BY OUTeUT. TG 
TEMPORARY DISK SPACE FOR USE IN SUBSEQUENT CALCULATIONS 
AND/OR LISTING. UNDER *MTS* TEMPORARY FILES NEED NOT 
BE YCREATED*® BUT THEIR CONTENTS ARE DESTROYED AFTER 

Ss. GNC Es 


TO RUN THIS PROGRAM FROM THE TERMINAL», THE FOLLOWING 
DEVICES MUST BE SPECIFIED IN THE RUN COMMAND: 


UNIT C MAGNETIC TAPE DEVICE (EG. *T*) 
4 OUTPUT TEMPORARY FILE FOR FURTHER CALCULATIONS 
5 (DEFAULTS TO *SOURCE*) 
6 (DEFAULTS TO *SINK*) 
8 OUTPUT TEMPORARY FILE FOR LISTING OF RESULTS 
9 INPUT FILE OF X-COORDINATES»s FLAP. ETC. 
EXAMPLE: 


SRUN SCANREADOB O=*T* 4=-S1 8=-L1 9=FX61162 


TO RUN THIS PROGRAM *BATCH*¥s THESE ASSIGNMENTS MUST BE 
CHANGED TO: 


UNIT 5S DISK FILE OF CALIBRATION PARAMETERS 
8 *SINK* 


EXAMPLE : 
SRUN SCANREADOB O=*T* 4=-S1 S=CALFILE 8=*SINK* 9=FX51163 


AFTER SEVERAL RUNS OF SCANREADOBs OUTPUT DATA ON UNIT 4 
CAN BE COLLECTED INTO ONE FILE BY USING THE INSTRUCTION: 


SCOPY -S1+-S2t+..-+-S9 TO SCANDATA 


THEN FILE SCANDATA BECOMES INPUT TC SCANPLOT ANDO 
SCANLIFT. HOWEVER» THE LAST LINE IN SCANDATA MUST BE® 


END 
TO OBTAIN A COMPLETE LISTING OF RESULTS: 
SCOPY =—L14-Let...t-LS TO SCANLIST 


$LIST SCANLIST 


INTEGER*¥2 IN(40) »sLENeL ABC 80) 
INTEGER NAS(50)+e+NSCAN( 50) »sWORK( 200) sNAME(5) »>NPST (S50) »«NTP( 50) 
DIMENSION Q(5C) »SCAN( 50) »CP(50).X(5C) sPST(SO) »TP(S5S0) 


READ IN THE LABEL BLGCK FROM TAPE 
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CALL READ(CLABsLENs0eLNR,0,6100) 
CALL CONVAE(LABeLEN) 
WRITE( 6.6000) 

WRITE (626001) (LAB( I) »I=1+eLEN) 
WRITE( 6» 6002) 

READ(S+5000) ALPHA 
WRITE(6s6003) ALPHA 


READ IN THE FIRST DATA BLOCK 


CALL READCINsLENs OsLNR,05E101) 
LEN2=LEN/Z2 


INITIALIZE SUBSCRIPTS FOR DATA SEARCH 


Ka 

DOlS) L=26LEN2] 
WORK (K )=IN(L) 
K=K+1 

CONT INUE 


READ IN THE REMAINING DATA BLOCKS 


CALL READCINeLEN20sLNR20,E102) 
LEN2=LEN/2 

DD Se LH=1S LENS 

WORK (K)=IN(L) 

K=K+1 

CONT INUE 

GO TO 4 


LEN2=LEN/2 

D016, L=145LEN2 
WORK (K )=IN(L) 

K=K4¢] 

IF(K.GT.200) GO TO 7 
CONTINUE 


READ INPUT DATA FROM DOISK FILE 


READ(9,9000) NPTSsFLAPs(NAME( I) » 1=1.5) 
READ(9+9001) NMUsNML oo NFU oNFL 
READ(9.9002) (X(1I)+I=1sNPTS) 


ASSIGN COUNT VECTORS 


Ja 

CUS ll saNers 
NAS(I)=wORK(J) 
NPSTCI )=WORK(J+1) 
NSCANC 1) =WORK(J+2) 
NTP( I) =wORK(J+3) 
J=JIt4 

CONTINUE 
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CUTPUT *tHE<«GOUNT VECTORS 18 DISK FOR LISTING 


WRITE(8s8000) (NAMECI)»,I=1,5) 

WRITE( 8s 8001) ALPHAsFLAP 

WRITE( 8, 8002) 

WRITE( 8s 8003) (C(CIseNAS(I) »sNSCANCT)SNPSTC(IT) sNTPCI)) »IT=1eNPTS) 


READ IN CALIBRATION CONSTANTS 


WRITE( 646007) 
READ(S5+5001) QSL»QINT 
WRITE(6,6008) QSL»QINT 
WRITE( 6, 6009) 

READ(5 +5001) SStL»sSINT 
WRITE( 626010) SSt»SINT 
WRITE( 65,6011 ) 

READ(5 .5001) PSL»PINT 
WRITE(6,6012) PSL»PINT 
WRITE( 6, 6013) 

READ(5 5001) TSL»TINT 
WRITE( 636014) TSLeTINT 


CONVERT RAW COUNTS TO PHYSICAL UNITS 


DO 9 I=leNPTS 
Q(1)=QSL*NAS(1I)#QINT 

PST(C 1) =PSL¥*¥NPST(I)4#PINT 
TPC I)J=TSL¥NTPCI)4TINT 
SCAN(I )=SSL*XNSCAN( I) #SINT 
CP( I )=1.C+#SCAN(I)7Q(1) 
CONT INUE 


COPY VECTORS TG DISK FOR LATER) LISTING 


WRITE( &s 8004 ) 
WRITE( 8 8005) (€15Q(1I)+SCAN(I),CPC(I)»-PSTC I)» TPC(I)) sIT=1eNPTS) 
WRITE( 8s 8006) 


QUTPUT RESULTS TO DISK FOR FURTHER CALCULATIONS 


WRITE( 4.4000) (NAME(I)>,I=1-5) 

WRITE( 4, 4001) ALPHAsFLAPs NPTS eNMU oNML oe NFUs NFL 
WRITE( 454002) (X(1I)+I=1+NPTS) 

WRITE (4.4003) (CPCI) + I=1+NPTS) 

WRITE 65 6004) 

GO TO 9399 


FORMAT (SAG) 
FORMAT (2F4.02513) 
FORMAT (1 6F 5.4) 
FORMAT (1 OF 8. 3) 
FORMAT (F 4.0) 
FORMAT (2F8.3) 
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6000 FORMAT(/S*FILE LABEL? *) 


6001 FORMAT (/740A2) 

6002 FORMAT(/*ENTER ALPHA: *) 

6003 FORMAT(/S* CONTINUING WITH ALPHA=",.F5.0) 

6004 FORMAT (/*DATA COPIED TO DISK FOR FURTHER CALCULATIONS *) 

69005 FORMAT (/*END-OF=-FILE ENCOUNTERED IN LABEL BLOCK") 

6C06 FORMAT(/*END-OF-FILE ENCOUNTERED IN FIRST OATA BLOCK® ) 

6007 FORMAT(/Z*ENTER Q@ CALIBRATION SLOPE AND INTERCEPT: ®*) 

6008 FGRMAT(* CALIBRATION IS: =",F8.59**NASt® F522) 

6009 FORMAT(/Z*ENTER SCAN CALIBRATION SLOPE AND INTERCEPT: ® ) 

6C10 FORMAT(*CALIBRATION IS: SCAN=",F8.50 **NSCAN+t* 5F 522) 

6011 FORMAT(/*ENTER PST CALIBRATION SLOPE AND INTERCEPT: *) 

6012 FORMAT (* CALIBRATION IS: PST="sF8.4e"*NPSTt+* F621) 

6013 FORMAT(/*ENTER TP CALIBRATION SLOPE ANO INTERCEPT? ® ) 

6014 FORMAT(*'* CALIBRATION IS: TP="* sF8.45"*NTP+* »F5.1) 

800C FORMAT (//* DATA FOR TEST OF" 92Xe5A4) 

8001 FORMAT(/S*TEST CONDITIONS: ALPHA="sFE.0e*% FLAP="»F6.0) 

8002 FORMAT (/*RAW COUNTS#/T351HI sT8s3HNAS 3718s S5SHNSCANo T2998 
14HNPST »6T40.3HNTP) 

8003 FORMAT (13,4110) 

8004 FORMAT (/*PHYSICAL QUANTITIES**/4HHOLE 2T10s1HQsT20+s4HSCAN » 
1730+ 2HCP o6T38%3HPST »T48s1HT) 

BOOS FORMAT (1 39 3X% sF 70 39 3X 9F le 30 3X oF Velo 3XeF e290 4XeF4el ) 

8C06 FORMAT (47 ** QUANTITIES NOT REQUIRED ARE MASKED ***** x ®) 

900C FORMAT (135F4.0 »5A4) 

9001 FORMAT(414) 

9002 FORMAT (41 6F 5.4) 


(c SREGLALS EA LTS 


100 WRITE(6,6005) 
GO TO 9399 
101 WRITE 66,6006) 


999 STOP 
END 
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PROGRAM NAME: SC ANPLOT 


OBJECT NAME: SCANPLOTOB 


THIS PROGRAM READS DATA WRITTEN TI DISK BY SCANREAD AND 
PEOTS The PRESSURE DISTRIBUTION FOR EACH SE OF DATA. 
NOTES TERMINATION OCCURS WHEN * END*® IS ENCOUNTERED IN 
THE DATA FILE ON UNIT 5 (THE INPUT DATA) IN DATA LINE 
ONE. I/O UNITS SPECIFIED IN THE RUN COMMAND ARE: 


UNIT 5S FILE OF INPUT DATA 
6 (DEFAULTS TO *SINK*) GUTPUT LIST 
9 TEMPORARY DISK FILE TO BE USED WITH *CALCOM2aQ 


THE PLOTTER SUBROUTINE LIBRARY MUST SE USED. 
A TYPICAL RUN CCMMAND WOULD BE: 
SRUN SCANPLOTOBtXPLOTLIB S=SCANDATA 9=-TEMP 


IF RC=C ON RETURN FROM SCANPLOTOBs THE PLOTTING IS 
GOMPEETED BY: 


SRUN *CALCOMPQ SCARDS=-TEMP 


DIMENSION WORK(2C48) »X(50).CP(50) 

INTEGER END/Z*® END*Z,INZ EP Of * 4 sALPHAsFLAPsNAME(5) 
CALL PLOTS( WORK »81°2) 

CALL PLOT(2.022.0,-3) 

COUNT=1- 0 

FACT=0 .667 


SFACT OR * “SUBROUTINE, SCALES TEE PHYSICAL. SIZE OF THE 
RESULTING PRESSURE DISTRIBUTION PLOTS. FACT=1 WILL 
PRODUCE A PLOT 1C INCHES BY 6 INCHESs THE MAXIMUM 
ALLOWABLE SIZE FOR THIS PROGRAM, 


CALL FACTOR(FACT) 
READ INPUT DATA 


READ(S5»S5000) (NAMECI)+s1=155) 
IF(NAME(1).=QeEND) GO TO 999 

READ(S 45001) ALPHAsFLAPsNPTS »+NMU 9 NML » NFU» NEL 
WRITE( 6s 6CC0) ALPHA 

READ(S5 .50C2) (X(1)sIT=1s+NPTS) 

READ(54£903) (CPCI) »2I=1sNPTS) 


PLOTTING OF AXES AND LABELS 


CRUL LA KP SO6620s 0 6 C's 2ZHCP = 2 56 69'20605=5.001.C810.0) 
CALL AXIS(05.0 20.0 91HX9-1510266290 .0+-0.020.1510.C) 
CALL PLOT(5.020.09-3) 

DOP2 l=", 5 

Y=FLOAT(I)Z1.2 
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CALL SYMBOL(=-5.CoY2C.250 »sNAME( I) »90.0 24) 
CONT INUE 


CALL SYMBOL (=-4 ..501.050.125s EHALPHA= 290.036) 
CALL SYMBOL (-4.592.02C.125sALPHA290.C 94) 
CALL SYMBOL (=-4 6.0 o1.05C.125s SHFLAP=s 990.035) 
CALL SYMBOL (-4.052.020.125sFLAP290.6 54) 


IF(FLAP.ZEQ.IN) GO TO 20 
PLOT SUBROUTINE FOR FLAP=0UT 


CALL FLOUT(X+«CPs+NMU sNML » NFU SNFL ) 
GO TO 40 


PLOT SUBROUTINE FOR FLAP=IN 
CALL FLIN(XsCP»NMUsNML ) 
SET UP PLOT ORIGIN FOR NEXT PLOT 


PF CGOUNT st 60.0) GO. T@n60 
CALL PLOT(-5.0912.0.-3) 
GO TO 7C 

CALL PLOT(7.02-12.0s-3) 

C OUNT= -C OUNT 

GOeTO7 4 


FORMAT (5A4) 

FORMAT (244.513) 

FORMAT (1 6F 5.4) 

FORMAT (1 OF 8. 3) 

FORMAT (/*PLOTTING FOR ALPHA=* 462XsA4) 


PROGRAM EXIT 


TRCCOUNT SGT. 0.0) “GO TC, 80 
GALE  PEOTC(.0s=12.093—3) 
CALL PLOT(0.C20.0+999) 
STOP 

END 


SUBROUTINE FLIN(X»CP »NMU »NML) 
THIS SUBROUTINE PLOTS FOR FLAP=IN CASE 


DIMENSION X(5C) »CP(50) »XGR(30),CPGR(30) 
UPPER SURFACE 


Dowd 1=1 .NMU 
XGRCI)=xX(1) 
CPGR(I )=cCP(T) 
CONT INUE 
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CALL DIST(XGReCPGR oNMU 91 ) 
LOWER SURFACE 


DO 2 1=1 sNML 

J=NMU+I 

XGR( I) =X (J) 

CPGR(I )=CP(J) 

CONTINUE 

CALL DIST(CXGRsCPGReNML 2) 


RETURN 
END 
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SUBROUTINE FLOUT(X»CPsNMU»sNML »NFUSNFL ) 
THIS SUBROUTINE IS FOR FLAP=0UT “CASE 
DIMENSION X(50) sCP(50) +XGR(20)sCPGR(20) 


a) 


J 


MAIN UPPER SURFACE 


DO 1 I=1.NMU 

XGR(1)=x(1) 

CPGR(I )=CPC(1) 

CONT INUE 

CALL DISTCXGReCPGReNMUs1l) 


MAIN LOWER SURFACE 


OQ 2 I=1 »NML 

J=NMU+ I 

XGR(I)=xX(J) 

CPGR(I )=CP(J) 

CONT INUE 

CALL DIST(XGRsCPGRsNML 52) 


FLAF UPPER SURFACE 
DOS L=1 «NFU 
J=NMU4NML4I 
XGR(1I)=x(J) 
CPGR(I )=cCP(J) 
CONTINUE 
CALL DIST(XGReCPGRsNFUs1) 


FLAP LOWER SURFACE 
OO) 4° 1=1.sNFL 
J=NMU+4NML+4NF UFI 
XGRCI)=X (J) 
CPGR(I )=CP(J) 
CONTINUE 
CALL DIST(XGRsCPGRsNFL 92) 


RETURN 
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END 


SUBROUTINE DIST(XsCPsNsK) 
DIMENSION X(3C) ,CP(30C) 
X(N#1)=0.06 

X(N+2)=090.1 

CP(N+41 )=0.C 

CP(N+#2 )=1.0 

CALL FLINE (CPsXo—-No1515K) 
RETURN 

END 
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PROGRAM NAME: SCANLIFT 
OBJECT NAME: SCANLIFTOB 


THIS PROGRAM USES THE DATA PRODUCED BY SCANREAD TO 
CALCULATE THE LIFT AND MOMENT COEFFICIENTS. Wile 


INPUT DATA IS THE SAME AS THAT FOR SCANPLOT. THE 
PROGRAM MAY BE RUN THROUGH *BATCH* OR THROUGH THE 
TERMINAL. 


THE PROGRAM REQUIRES SUBROUTINES FROM *SSPMATHLIB. 


A TYPICAL RUN COMMAND WOULD BE: 


SRUN SCANLIFTOB+*SSPMATHLIB S=SCANDATA 


INTEGER END/* END*/,IN/* 0."/»sNAME(S) 
DIMENSION X(50) »CP(50) 


READ DATA DEPOSITED BY SCANREAD 


READ(5.5000) (NAME(1),I=1 25) 
IF(NAME(1).EQ.END) GO TO 999 

READ(5.5001) ALPHAsFLAPsNPTS+NMU oNML » NFUSNFL 
READ(S5+5002) (X(1I)»I=1sNPTS) 

READ(5.5003) (CP(I)»+I=1.NPTS) 

WRITE( 6.6000) (NAMECI) so I=125) 

WRITE( 626001) ALPHA»sFLAP 


IF(FLAP.NE.0.0) GO TO 2 
INTEGRATION SUBROUTINE FOR FLAP=IN CASE 


CALL FLIN(X+CP»NPTS »NMUsNML oCNo CMLE ) 
GOpg oes 


INTEGRATION SUBROUTINE FOR FLAP=OQUT CASE 
CALL FLOUT(XeCP»sNPTS »NMU sNML »NFU »NFL»> CNo CMLE) 
CALCULATE CL AND CM1/74 


CL=CN*COS( ALPHA*3.21416/180.) 
CMQ=CMLE+CL/4. 


OUTPUT THE RESULTS 
WRITE( 656002) CNeCMLE 


WRITE( 626003) CL»eCMQ 
GO Tor. 


S000 FORMAT (SA4) 
5001 FORMAT (2F4.02513) 
5002 FORMAT (1 6F 5.4) 
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FORMAT (1 OF 8.3) 
FORMAT(S/* TEST OF * 32X%55A4) 
FORMAT (* ALPHA=*,F6.092X%o *FLAP=! »F 6.0) 
FORMAT (*CN=* »F7.2/°CMLE=* oF7. 3) 
FORMAT (* CL=* 6F7.2/°"CM1/4="* oF 7.3) 


STOP 
END 


SUBROUTINE FLIN(XsCPsNPTS »eNMU eNML oCNoCMLE ) 
DIMENSION X(50).CP(50) sX1(51) sCPI(51) »WORK( 100) +ARG(8)> 


1VALC8) »XUC25) »CPU(C 25) »XL(25) »CPL(25)5,Z(51) 


H=0.02 


SET UP AUXILTARY VECTORS 


OOF 1 VreH1 551 
XI (1)=(1-1)*H 
DO 2 I=1 eNMU 
XUCTI)=XC1) 
CPU(I)=CP(T) 
DO 3 I=1 .NML 
J=NMU+I 
XL(C I )=X( J) 
CPL(1) =CP(J) 


INTERPOLATE CP VECTORS AND INTEGRATE FOR UPPER SURFACE 


LOOP=4 

GGa4 eh=1 «51 

CALL ATSG(XI (1) sXUsCPUsWORK »NMU ol sARGs VAL 58) 
CALL ALI(XI(1I)+ARGeVALsCPI( 1) »8sEPSsIER) 
IFC IER.EQ.C) GO TO 4 

WRITE( 66000) LOOPs,I1,IER 

CONT INUE 

CALL QSF(HeCPI,Z251) 

SUML UP =Z (51) 

O00 5S I=1.51 

CPI(I)=CPICI) *XICI) 

CALL QSF (HeCPI+Z251) 

SUMMUP =Z(51) 


DO LIKEWISE FOR LOWER SURFACE 


LOOP=6 

Pore t=) 551 

CALL ATSGCXI (I) »XL+CPL»WORKeNML»s1sARGoVAL 28) 
CALL ALI (XI(1)+ARGeVALsCPI( I) »8,EPSsIER) 

IF( IER.EQeC) GO TO 6 

WRITE( 65,6000) LOOPsI,IER 

CONT INUE 
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SUMLLO=Z (51) 

O17 Wein Ss 

CPIN I=CRI1) *¥x1'C1) 

CALL QSF(HeCPI.Z.51) 

SUMMLO=Z(51) 


COcCvVve CTVIRESULTS 


CN=SUMLL O-SUML UP 
CMLE=S UMML O-SUMMUP 


FORMAT(* IN FLIN LOOP *'sI5s* AT POINT* 145% TER=".14) 


RETURN 
END 


SUBROUTINE FLOUT(X+CPsNPTS»NMUsNML»NFUsNFL 2» CNoCMLE) 
DIMENSION X(€(5C) »CP(S0)+X1(40) sCPI1(40) »WORK(10C)+ARG(8)> 
1VAL(8) »XMU(30)+XML(30).XFUC20).XFL( 20) »CPMU( 20) »CPML( 30)» 
2CPFUC20) »CPFL(20).Z(50) 

EPS=1.0E-2 
H=0.02 


MAIN SURFACE 


OORiat=1..32 
XI (I )=(C1I-1)*H 
DO 2 I=1 +NMU 
XMUC IT) =X (1) 
CPMU(CT )=CP(I) 
DO 3 I=1 +NML 
J=NMU4 I 
XML (1) =X (J) 
CPML (I )=CP (J) 


INTERPOLATE CP VECTORS AND INTEGRATE 


LOOP=4 

DeRGel=1.32 

CALL ATSG(XI (1) »XMUsCPMU »WORK »NMUo1+ARGo VAL 58) 
CALL ALI(XI(CI)+sARGesVAL»CPI( I) »s8sEPS,IER) 
IFC IER.EQ.0C) GO TO 4 

WRITE( 6,6000) LCOPs,I>,IER 

CONTINUE 

CALL QSF(H»CPI+sZ2+32) 

SUMMUP =Z (32) 

BoOeSsy b=) 432 

CPI(1)=CPI(1I)*xXI(1I) 

CALL QSF (HeCPI 92932) 

ADDMUP=Z (32) 
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LC OP=6 

DO 6 TP=1532 

CALL ATSG(XI (1) +XML+CPML «WORK »eNML ol eARGe VAL » 8) 
CALL ALI(XICI)»sARG»VALsCPI( 1) +8sEPSs,IER) 
IFCIER.EQ.C)} GO TO”’6 

WRITE( 626000) LCOP.I»sIER 

CONTINUE 

CALL QSF (HsCPI+Z.32) 

SUMMLO =Z (32) 

OO" 7 “f=1.32 

CPI(CI)=CPICI)*XI(I) 

CALL QSF(HsCPIs 2232) 

ADDMLO=Z (32) 


FLAP 


EPS=0.1C 

DO 8 I=1.NFU 
J=NMU+4NML4I1 
XFUCT)=x(J) 
CPFU(I )=CP(J) 
00 9 I=1.,NFL 
J=NMU4NML4+NFU4I 
XFLCT)=x(J) 
CPFL(I )=CP (JU) 
DO" eho. Ta 14177 

XI (1 )= (1-1 )*Ht0.65 


INTERPOLATE CP VECTORS AND INTEGRATE 


LoOOP=11 

Deri a f= tis) 7, 

CALL ATSG(XI (1) +s XFUs CPFUs WORK sNFUel »ARGo VAL s NFU) 
CALL ALI (XI(1I) sARGsVALsCPI(I) sNNFUsEPSsIER) 
IF( IER.EQ.0) GO TO 11 

WRITE( 66000) LOOPsI>,IER 

CONTINUE 

CALL QSF(HsCPI+Z217) 

SUMF UP =Z (17) 

bo Je i=1, 17 

CPICI)=CPI(I)*XI(I) 

CALL QSF(HseCPI+Z:17) 

ADDF UP=Z(17) 


LOOP=13 
OO) is T=Ts17 


CALL ATSG( XI (1) +XFL»CPFL «WORK »sNFL 51 »>ARGs VAL »NFL) 
CAL ALI(XI(1)+ARGsVALeCPI( 1) eNFL»EPS»IER) 

LEP CIER .EG.0).GO 10 1s 

WRITE( 626000) LOOPsI,IER 

CONT INUE 

CALL QSF(H»eCPI+Z+17) 

SUMFLO=Z (17) 


Soca l= 1517 
CPI(C1I)=CPI (I) *XI(1) 


CALL QSF(HeCPI+Ze17) 
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ADDOFLO=Z (17) 
COLLECT RESULTS 


CN=SUMML O- SUMMUP+ SUMFLO-SUMFUP 
CMLE=ADD ML O- ADOMUP+ ADDFLO-ADDFUP 


FORMAT C* IN FLOUT LOOP".14s2Xs*AT POINT'51IS,? 


RETURN 
END 
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PROGRAM NAME: TRAVREAD j 


OBJECT NAME: TRAVREADOB 


THIS PROGRAM READS MAGNETIC TAPE RECORDS FOR TRAV MODE 
OPERATION AND REDUCES THE RAW COUNTS TO PHYSICAL UNITS. 
THIS IS FOLLOWED BY OUTPUT TO TEMPORARY DISK SPACE FOR 
USE IN SUBSEQUENT CALCULATIONS. THE PROGRAM MAY BE RUN 
FROM THE TERMINAL. HOWEVER, TO MINIMIZE I/70 THROUGH THE 
TERMINAL», RESULTS ARE COPIED TO TEMPORARY DISK SPACE FOR 
LATER LISTING ON THE LINE PRINTER. THE '"RUN* COMMAND 
MUST CONTAIN REFERENCE TO THE FOLLOWING UNITS: 

UNIT 0 MAGNETIC TAPE (EG *Tx) 
OUTPUT DISK FILE FOR SUBSEQUENT CALCULATIONS 
(DEFAULTS TO *SOURCE*) 
(DEFAULTS TO *SINK*) 
OUTPUT DISK FILE FOR LISTING OF RESULTS 
INPUT DATA FILE (SAME AS FOR SCANREAD) 


L$ 


ooonuw 


A TYPICAL RUN COMMANO FROM THE TERMINAL WOULD BE: 

SRUN TRAVREADOB O=*T* 4=-T1 8=-L1 S=INDATA 

A TYPICAL RUN COMMAND THROUGH *BATCH* WOULD BE: 

SRUN TRAVREADOB O0O=*T* 4=-T1 S=CALFILE 8=*SINK*® 9=FXK61163 


AFTER SEVERAL RUNS OF TRAVREADOBs DATA OUTPUT ON UNIT 4 
MAY BE COLLECTED INTO ONE OISK FILE FOR TRAVORAG USING: 


SCOPY -T1t-T2tee.et-T9 TO TRAVDATA 


THENs TRAVDATA BECOMES INTPUT DATA FOR TRAVDRAG. AE 
LAST LINE IN DRAGCATA MUST BE * END* TO TERMINATE THE 
EXECUTION OF TRAVCRAG PROPERLY. 


A CCMPEETE CISTING, OF RESULTSZ MAY SE OBTAINED BY 


S$COPY -L14-L2teeet-L9 TO TRAVLIST 
$LIST TRAVLIST 


INTEGER*2 IN( 40) eLENeLAB( 80) 

INTEGER NAS(50) seNWAKE(5C ) »sNPOS( 50) »WORK( 200) » NAME(5) » 
1NPST(5C) »sNTP(5O) 

DIMENSION Q(SO) «WAKE( 50) »POS(50) »PST(50)+TP(S5C) 


READ THE LABEL BLOCK 


CALL READCLABsLENs 0 oLNR2046100) 
CALL CONVAE(LAB>eLEN) 
WRITE( 6% 6000) 

WRITE(636001) (LABCI) »sI=1+LEN) 
WRITE( Es 6002) 
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READ(5+5000) ALPHA 
WRITE( 6s 6003) ALPHA 


READ IN THE FIRST DATA BLOCK AND PERFORM DATA SEARCH 


CALL READCINeLENeOsLNRoC.E101 ) 
LEN2=LEN/Z2 

K=1 

DO_ 3, ib=2.LENe 

WORK (K J=IN(L) 

K=K+1 

CONTINUE 


READ IN THE REMAINING DATA BLOCKS 


CALL READ(CINeseLENsO0sLNRO 566) 
LEN2=LEN/2 

DO 5 L=1,LEN2 

WORK (K )=IN(L) 

K=K+4+1 

IF(K.GT.~20C) GO TO 8 

CONT INUE 

GO TO 4 


LEN2=LEN/Z2 

DO 7 L=1 .LEN2 
WORK (K )=IN(L) 

K=K+1 

IF(K.GT.200) GO TO 8 
CONT INUE 


READ INPUT DATA FROM DISK 
READ(9+29000) NPTSeFLAPs(NAME(C I)sI1=1.5) 


ASSIGN COUNTS TO WORKING VECTORS 


K=K/5 

jai 

DO_ 9. V=1 sk 
NASC(1I)=WORK( J) 
NPST (I )=WORK(J+1) 
NWAKE( 1) =WORK(J+42) 
NPOS (1 )=WORK(J+3) 
NTP CI) =WORK(J5+4 ) 
J=JI45 

CONT INUE 


OUTPUT COUNT VECTORS TO OISK 


WRITE( 8s 8000) (NAMECI) -T=155) 

WRITE( 8+ 8001) ALPHA»sFLAP 

WRITE( 8, 80C2) 

WRITE 8s 8003) (CI sNAS(I) sNWAKECI) eNPOSC IT) eNPSTC(I) NTPC I) ) »IT=15K) 
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READ IN CALIBRATION CONSTANTS 


WRITE (6s 6007) 

READ(5 +5001) QStL»QINT 
WRITE(E€+6008) QSL+QGINT 
WRITE(C 6, 6009) 

READ(5 2.5001) WSL».WINT 
WRITE( 6.6016) WSL»WINT 
WRITE( 646011), 

READ(5 +5001) PSL»PINT 
WRITE( 656012) PSL»PINT 
WRITE( 6, 6013) 

READ(5 +5001) PSTSL»PSTINT 
WRITE( 6s 6014) PSTSL»sPSTINT 
WRITE( 646015) 
READ(S+.5001) TSt»TINT 
WRITE(6,6016) TSL»sTINT 


CONVERT COUNTS TO PHYSICAL UNITS 


DO 10 I=1+K 
Q(1I1)=QSL*NAS(I)4QINT 

WAKE CT )=WSL*NWAKE (I) ¢WINT 
POS(I)=PSL*NPOS(I)4#PINT 
PST(1)=PSTSL¥NPST(I)+PSTINT 
TPC I )=TSL*¥NTP(T)+TINT 
CONTINUE 


COPY RESULTS TO DISK FOR LATER LISTING 


WRITE( 8s 8004) 
WRITE( 8s 8005) (C1eQ(1) sWAKECI)sPOS(I)+sPST(I)+TPCI))»I1=1+K) 


WRITE( 8s 8006) 
QOUTPUT RESULTS FOR FURTHER CALCULATIONS 


WRITE( 424000) (NAMEC(I)+.1I=155) 
WRITE(4s4001) ALPHAsFLAPskK 
WRITE( 4.4002) (WAKECI) »I=1+K) 
WRITE( 4, 4003) (POS(I)+IT=1+K) 
WRITE(4,4002) (Q( 1) eIT=15K) 
WRITE( 656004) 

GO TO 999 


FORMAT (5A4) 

FORMAT (2F4.0 »13) 
FORMAT (1 OF 864) 
FORMAT (1 6F5.3) 
FORMAT (F 4.0) 

FORMAT (2F8.3) 

FORMAT (/*FILE LABEL: *) 
FORMAT (/40A2) 

FORMAT (/*ENTER ALPHA®*) 


a ee nr . ce a ts - . 
- ‘¥ 4" ‘\ 


a a : 


ae 7 > Gan 
. ae a er MATene> ova Z 


-_ 7 Ger oe ; 
: THrOe is as 
ion TH1D6 kanes atte 


ratw, Jaw (vor, = aKa 
Wty jew Lo1caw8 a7 ; 

: 1 115000 
THISyJe9 (1008, 2% 
nr. f20 (sass STOR 
(erodes )aTIAW 
TUIFea, ueT2e (sees. 2PQnea 
Tatt2heo2Ted (ei.ce .o eT iow 
talon. pet taw 
fHtt+det (1908<2}0aRR 
THITeLISl LSEte4a.4Id3TiVNe 


>TTAU JAD leyHde UT ema? rAaveag?> 


nate? cing 
TrAIVe CT TEA Me eet) 
*Hine Di TSanwhe Jewel 1) 
TH 4a! C3 2c4uMe.s* a1 209 — 
Tv (29407 OT 2 ahe 2 feR=41 F248 
rviTed petueyeret: tat 
Ay SUM THD 
‘jy 
DALTZEL. RATAY ACF HATE OF eT vac9 


: 1a0G8 ASTIAW 
4s Pele (11 79S9Te bP) Teel 1 1209111 esAws Toni) (2508 6a STIAW 
"= (2008.8 1 STinW 7 


2VO1 TAUUDUAS SSRTAUA BON <Tguese rusTyO 


=te(t)aMauM) (o0ne ss dSviRe 

Ws nd he 24 (1000.8 197TH 
(As i= te tT DaaAW) (sdge.a)sTiAW 
(Hef=lelt 2204) PCOOA,s VETIAW 
Cis Pte bE [3698,0j)5TI9W 
a Bee 


B19 
74 
6003 FORMAT (/Z*CONTINUING WITH ALPHA=* ,F6.0) 
6004 FORMAT(S*DATA COPIED TO DISK FOR FURTHER CALCULATIONS.*) 
6CC0S FORMAT (/*END-OF-FILE ENCOUNTERED IN LABEL BLOCK.*) 
6006 FORMAT (/*END-OF-FILE ENCOUNTERED IN FIRST DATA BLOCK.") 
6607 FORMAT(/*ENTER Q@ CALIBRATION SLOPE AND INTERCEPT: *) 
6008 FORMAT(*CALIBRATION IS: =*,F8.5e" XNASt* F522) 
6COS FORMAT(/*ENTER WAKE CALIBRATION SLOPE AND INTERCEPT: *) 
6010 FORMAT (* CALIBRATION IS: W=",s9F8.5e" *NWAKE+',FS.2) 
6011 FORMAT(/*ENTER POS CALIBRATION SLOPE AND INTERCEPT: *) 
6012 FORMAT(* CALIBRATION IS: POS=",F8.55"**NPOSt'5F5.2) 
6013 FORMAT(/S*ENTER PST CALIBRATICN SLOPE AND INTERCEPT: *) 
6014 FORMAT (* CALIBRATION IS? PST=",F8.45"*NPST+"*sF6.1) 
6C15 FORMAT(/*ENTER TP CALIBRATION SLOPE AND INTERCEPT? ®*) 
6016 FORMAT(* CALIBRATION IS: TP=*,F7.25"*NTPtH! s F662) 
8000 FORMAT(//"* DATA FOR TEST OF * 42X%s5A4) 
8001 FORMAT(/S*TEST CONDITICNS: ALPHA=",FE6O.0s* FLAP="sF6.C) 
8002 FORMAT (4 *RAW COUNTS*/T391HI 9 T10s3HNASoT19 sSHNWAKE oT 3C 
14HNPOS +740 s4HNPST TSC » 3HNTP) 
8003 FORMAT (13,5110) 
8CO04 FORMAT (/* PHYSICAL QUANTITIES**//T 351HI+T1091HQsT20s 4HWAKE 9 
1730+ 3SHPOS»s T40¢e 3HPST »T505 2HTP) 

8COS FORMAT (145F10.292F10.394X oF 64.194XsF 60.1) 
B8CCE FORMAT (/** QUANTITIES NOT REQUIRED ARE MASKED ****x®) 
9000 FORMAT (1 35F4.0+5A4) 


cS SPECI AE EXITS 


100 WRITE(6, 6005) 
Gor To 999 


101 WRITE Es 6006) 
Cc 
G EXIT 
S 
999 STOP 
END 
Cc 
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1D 
PROGRAM NAME: TRAVDRAG 


OBJECT NAME: TRAVDRAGOB 
THIS PROGRAM CALCULATES THE DRAG COEFFICIENT FROM THE 
WAKE DATA PRODUCED BY TRAVREAD. IT MAY B= RUN THRIUGH 
THE TERMINAL OR THROUGH *BATCK*. 


IT REQUIRES SEVERAL SUBROUTINES FROM *SSPMATHLIB. A 
TYPICAL RUN COMMAND WOULD BE: 


SRUN TRAVDRAGOBt+*SSPMATHLIB S=TRAVDATA 
THE PROGRAM TERMINATES ON READING * END®* IN THE FIRST 


LINE OF A DATA SET. 


INTEGER END/*® END*/ +NAME( 5) 
DIMENSION WAKE(S0O)+Q(50) »sPO0S( 50) s0(50)5Y(106)sP(10C) >» 


1ARG(8) eVAL(8)+Z(100) »POSF (S50) 


EPS=1.CE—2 
H=0.005 


SET UP AUXILIARY POSITION VECTOR 


DO 1 IT=1+100 
Y¥(1)=C 1-1) *H 


READ DATA DEPOSITED BY TRAVREAD 


READ(5.5000) (NAMECT) sI=155) 
IF(NAME(1).EQeEND) GO TO 999 
READ (55001) ALPHAsFLAPsK 
READ(5559202) (WAKE( I) + I=1 5K) 
READ(5.5003) (POS(I).I=1>3K) 
READ(5.5004) (Q(1I)s1I=1+K) 


CALCULATE AUXILIARY VECTORS 
DO 3 I=1 6K 


D(I)=WAKECI)/ZQ(1) 
POSF (I )=POS(1I)/12. 


INTERPOLATE VALUES AT C.00S FOOT INTERVALS AND INTEGRATE 


Oe aarti t lOO 

CALL ATSGC(Y(I)sPOSFsD0sWORK Ke 1le2ARGs VAL» 8) 
CALL ALI(CY(1)+*ARGoVAL »P(1)s8sEPSs IER) 
IF(IER.EQ.0) GO TO 4 

WRITE(€, 6000) IsIER 

CONT INUE 

CALL QSF (HsP2Zs100) 


CD=Z(100) 


OUTPUT THE RESULTS 
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WRITE( EG, 6001) (NAMECT) e1IT=1-5) 
WRITE( 6.6002) ALFHAsFLAP 
WRITE( 6.6003) CD 

GO.i1 0 2 


FORMAT (SA4) 

FORMAT (2F4.0 213) 

FORMAT (1 0F 8.4) 

FORMAT (1 6F 5.3) 

FORMAT (1 OF 8.4) 

FORMATC*AT POINT*® »IS5* IER=",I5) 
FORMAT(/S/* TEST OF * 42X95A4) 
FORMAT (* ALPHA="4F6.0%" FLAP=" .F6.0) 
FORMAT (* DRAG COEFFICIENT=!,F9.5//) 
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APPENDIX C 


COMPUTER OUTPUTS 


Cl 
84 
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STUSTUN. AT ao 


SCON *T* POSN *14*% 
SR SCANREADOB O=*T* 4=-A 5=ZEROS 8=*SINK* 9=PXFZERO 
VTe2 7.43 


FILE LABEL: 


ENTER ALPHA: 


CONTINOING WITH ALPHA= 0. 


DATA FOR TEST OF FX-61-163 SERIES IT 


TEST CONDITIONS: ALPHA= 0, FLAP= 0. 
RAW COUNTS 

oi. 2awg RAS __ NSCAN NPST NTP 

1 467 244 29 27 

2 467 285 29 oa 

3 469 341 29 27 

4 468 S0g 29 2 

5 467 396 29 28 

ae 468 ei 29 2 

41 466 4 30 29 28 

8 466 437 29 ely 

9 470 426 29 27 

10 468 HOC 29 27 

14 471 367 29 28 

moe oo. 3A M 29 per 27s 

ee 467 340 29 27 

14 467 30C 29 27 

15 469 2U6 29 28 

16 467 350 29 te | 

17 470 362 30 27 

y o18 eso 6G 29 27 

19 469 378 29 28 

20 468 387 29 28 

21 4u72 381 30 27 

22 469 375 29 pop 

23 466 346 29 oF 

24 069 2 357 29 28 

= aS 315 29 28 

26 470 305 29 28 


ENTER Q CALIBRATION SLOPE AND INTERCEPT: 
CALIBRATION I[S: y= -0.08240 *NAS+59, 88 


ENTER SCAN CALIBRATION SLOPE AND INTERCEPT: 
CALIBRATION IS: SCAN=-0. 20200*NSCAN+ 38.70 


ENTER PST CALIBRATION SLOPE AND INTERCEPT: 
GAL LORATION L590) (PS TH*teeeeCEeeECN PST Sees ee 


ENTER TP CALIBRATION SLOPE AND INTERCFPT: 
CALIBRATION IS: TP=** #2 HK NIT D+ i ok ok 


PHYSICAL QUANTITTES* 
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jie 90 10. 588 0.51 ROK tok 
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ipa WA healer £189 THT 
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eel a cree ® _ 
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eK take keke ak 


0.42 ee kk 
3 73a Pipe WWE esc Ve Gh. -0.42 Ra tok teak ok C3 
4 Zalines | 7 -36.444 -9.71 a OK tok ok we ee ak 26 
5 oa a ae =O 7,292 =o ek ke ak ocak ok te 
6 fe APE ie I) -46.746 -1.19 eke HI ok ake ok ok 
— re ere. -48. 160 ole Fp | HK RK KK a Kr 
8 Pa oth “49,574 -1,31 Hee ee hea fe ae HC ok 
< Lagos =U 7 32 -1.24 We aR HK a ak fe kK 
10 eluate leh -42.100 -0.97 Rok kok kc ak CK 
i Fae Oe 8s 8 3 esd ou -0.68 ek kak Oko OR 
2 Pee ese ie I =3A. 090 -0.49 FOR HOR OK eK a ak 
= ne eles hoe ee ee OU -0.40 fe ak kak a eK ek 
14 21.399 “24.906 cafes 8 KKK KES ek I ok 
15 2 Ave a0, -~10.992 0.48 eK ok a kk ok ee 3 
16 21399 =e PPE EOL: =0,550 tok kik ak we te cok 
17 Ziliemioe -34,42u -0.63 eRe kate ok ak ok ek ie ake 
18 ee MAP May = 3s ee =Du.67 Fe kK Fk aK Ok 
7 19 21.238 -37.656 =On er 3K ee ake 2k ak sek ok 
20 Ping  -—39.208 “0.85 kK tok Ok fo OK 
ou 20984 = 3S enue -0,,8.2 OK Ke KK 
a2 2 liemat Sse a) -0.74 Heke ea ake ok 2 eK 4 
oP) ee eats Be -31.192 -0.45 Fk HK Hk ak ec ak ak 
24 ZA eee =—32 eyed =O, 52 KKK KKK % 
20) 20.981 =25 5930 -0,419 OR RK KK Hk toe 
26 21 Cgow =—22.040 -0.08 Ae te ok a 4 XK 


* QUANTITIZS NOT REQUIRED ARE MASKED *****£*% 


DATA COPTED, TO DISK, FOR FURTHER CALCULATIONS 
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FILE LABEL: _ 


ENTER ALPHA 


CONTINUING WITH ALPHA= 0. 


DATA FOR TEST OF FX-61-163 SERIES TT 


C4 
87 


TEST CONDITIONS: ALPHA= 0. FLAP= ae 
RAW COUNTS 

Sees NAS NWAKE NPOS NPST NTP 
1 496 830 579 29 28 
2 4ug7 829 540 29 28 
3 4ug7 B29 506 29 27 
4 496 829 473 29 27 
5 ugg 829 44d 29 27 
6 it ee eee 414 29 Zi 
7 498 839 386 29 27 
8 499 B27 359 29 8] 
9 497 827 334 29 28 
10 497 824 308 29 Pe 
11 499 stale 284 29 27 
een a4. u96 64C ae 29 27 
13 4uo4 452 267 29 | 
14 4a9g 311 261 30 27 
15 496 112 254 29 oT 
16 498 32 ale 29 ae 
1a 497 31 239 29 28 
18 498 30 231 29 oe | 
19 498 100 224 29 27 
20 498 2U6 218 29 27 
oa 496 372 21> 29 27 
a2 498 519 206 29 27 
23 498 646 201 29 27 
ie ee doa es 729 193 29 28 
meee 500 789 186 29 Pa | 
26 498 813 179 29 oA | 
27 496 822 158 29 27 
28 495 82u 132 29 27 
29 500 823 105 29 27 
30 499 825 Tt 30 28 


ENTER Q CALIBRATION SLOPE AND INTERCEPT: 
CALLTBRATION 1S¢ y=-0.08240*NAS +59. 88 


ENTER WAKE CALIBRATION SLOPE AND INTERCEPT: 


CALIBRATION TS: W=-0.00500¥*NWAKE+ 4, 25 


ENTER POS CALIBRATION SLOPE AND INTERCEPTS 
CALIBRATION IS: POS=-0.01100*NPOS+ 6.16 


ENTER PST CALIBRATION SLOPE AND INTERCEPT: 


CALIBRATION IS: PST=*####e ex EN PST + *¥ #4 4 
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FNTER TP CALIBRATION SLOPE AND INTERCEPT: 
TD = 8 eK FOR RN DD 4 kK Ok 


CALIBRATION IS: 


C5 
PHYSICAL QUANTITIFS* 88 
T .@) WAKE POS PST TP 

1 On 0.700 lo af a8 Mee He Ke ak ak eK KO tk 
2 18,93 0, 105 OF 220 We te keke Ie ak tok Kk 
3 ER surges Des! nas 0.594 ek KK 3k ak tee ee 
4 ola pa de Oto 05957 OK Ka Fk ek eke ak 
ate 18.84 Oe 1 1.320 ek kk tek ek ake ke oe ok 
6 18.68 O, 415 1.639 Fe OOK ok fk a ke 
a 18.84 0.100 1.914 ok $k tek We FOR KK 
8 Ar eLo 0.1715 Pg 2 es | Fe KO a ok We ake ae kc ak 
9 eno 3 0.175 Fa Wal) ek eke a ak $e Fa ok 
10 18.93 0.130 eee ee te Re ako eo HK 
ee eee Ge 0.250 3.936 ek 4 OK OK 
12 79 Ou 1.050 331.57 ek tee eH ek ek 
13 gis TAG et 1290 Cli aesG | Ok koa a Sek ef 
14 13,75 2695 Bee 2K ae ake a ak a ke ok ak 
15 Mist aco: | 3.660 3.366 eK fea ok ak ee a 
16 18.84 4.090 3,443 fe eo ak oc ea ak aie oe 
et oe ee 10003) _4.095 3353) dated de ee eke a ak 
18 18.84 fe 414) 0 356.19 We ake ake ae He OK 
19 18.84 Chasse) 3.696 We fe ae 3K 9k $e ok oe ok ok 
20 18.84 3.020 3162 3 ek ak 3 ek ke ak 
21 13508 BEALS KE B5 B2n Hee te ak ak a Seok KK 
Ze 18.84 7,655 3.894 eo a Xe ke ak He kK 
123 oes 88 1.020 3.949 tok Ak 7K aK a 
24 Abe BR 0.605 4,037 kK ak ae a ak Fk ak ok ok ak 
Zo Eon Oo Gus 4.11714 Se He oR ok we ea ok ak 
26 18.84 0.185 4,191 BOR Fe kok $K tek ye ak 
0. ec 0.140 4,422 ee eo ae 3k ok dete kok 
28 1950S 0. 139 4,708 Fe fk a ok He teak a ok 
29 18.68 Oa 3 Sao sala sdiaidaaats eK OK 
30 ieee 0 125 5.313 Fok RA OK Fok I Fok 


* QUANTITIES NOT REQUIRED ARE MASKED ***** 


VALA COPIED TC DISK FOR FURTHER CALGULA TIONS. 
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